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In order to investigate the optical properties of rubrene we study the vibronic progression of 
the first absorption band (lowest ty — > 7r* transition). We analyze the dielectric function £ 2 of 
rubrene in solution and thin films using the displaced harmonic oscillator model and derive all 
relevant parameters of the vibronic progression. The findings are supplemented by density functional 
calculations using B3LYP hybrid functionals. Our theoretical results for the molecule in two different 
conformations, i.e. with a twisted or planar tetracene backbone, are in very good agreement with 
the experimental data obtained for rubrene in solution and thin films. Moreover, a simulation based 
on the monomer spectrum and the calculated transition energies of the two conformations indicates 
that the thin film spectrum of rubrene is dominated by the twisted isomer. 



I. INTRODUCTION 

Organic semiconductors have become very promising 
and actively investigated materials [TJ [21 E] . Due to their 
specific optical properties organic molecules are partic- 
ularly interesting for various device applications such as 
light emitting diodes, field effect transistors and solar 
cells jU [5j |6] . Different studies have demonstrated that 
the optical properties of molecular solids are strongly 
influenced by their structure. Crystalline materials as 
for example pentacene [7] and diindenoperylene [8 thin 
films exhibit complex and generally anisotropic absorp- 
tion spectra which are often influenced by characteristic 
charge transfer excitations. Amorphous materials with 
extended 7r-conjugated electron systems, as for example 
rubrene (Fig. [I]), show optical spectra which are seem- 
ingly simpler to understand [9l [TUJ [TT] . However, when 
considered in more detail the growth of amorphous ma- 
terials may present its own complications. For example, 
studies of rubrene thin films on silicon uncovered a sur- 
prising growth behavior: Using X-ray reflectivity rubrene 
thin films are found to be relatively smooth and amor- 
phous over a large thickness range [11 . Yet, an anoma- 
lous roughness evolution of the films is observed, i.e. a 
decrease of the film roughness at the beginning (below 
~20nm) and a subsequent increase of the roughness at 
a later stage. Recent near edge X-ray absorption fine 
structure (NEXAFS) results indicate the existence of dif- 
ferent rubrene isomers on Si02 and Au(lll) [10]. At the 
early stage of growth both species were identified, i.e. the 
rubrene with the twisted and planar tetracene backbone, 
followed by the planar isomer only. 

In the present study we focus on the optical proper- 
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FIG. 1: Structure of the rubrene molecule (C42H28) which 
consists of a tetracene backbone (TBB) and four attached 
phenyl rings. Due to the different orientation of the phenyl 
rings and TBB the 7r-conjugated systems of both parts are 
strongly decoupled. 



ties of rubrene molecules in solution and amorphous thin 
films. In particular, we investigate whether the different 
isomers can be identified in the optical spectra. For this 
purpose we calculate the optical response of the differ- 
ent isomers of rubrene with time-dependent density func- 
tional theory (TD-DFT). The deformation in the relaxed 
excited geometry of each isomer found with TD-DFT is 
projected onto its internal vibrations. This allows an 
assignment of the observed sub-bands in the absorption 
spectra to an effective vibrational mode which is an aver- 
age over the various internal vibrations elongated in the 
relaxed excited geometry. Furthermore, we use the calcu- 
lations to assess the contribution of the different rubrene 
isomers to the observed spectra. 
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II. EXPERIMENT 

The material used in our experiments was purchased 
from Acros and purified by gradient sublimation. This 
procedure yields small rubrene crystallites which are rel- 
atively inert against oxidation in air [9]. To avoid well- 
known thermochromic effects - rubrene changes its color 
towards the red at elevated temperatures [12] - we took 
all spectra at T = 25 °C. 

The absorption spectra of rubrene in acetone solution 
were recorded using a UV-Vis spectrophotometer (Varian 
Cary 50). With an overall reproducibility of the optical 
constants being on a 1% level we verified that for the 
chosen concentrations (1.46 x 10 -5 — 1.15 x 10 _4 mol/l) 
the absorbance data follow the Lambert-Beer law. The 
optical properties of rubrene thin films were measured in 
ultra-high vacuum using a spectroscopic rotating com- 
pensator ellipsometer (Wollam M-2000) with a broad 
band 75 W Xe lamp (250 - 1000 nm) and CCD based de- 
tection system with a spectral resolution of about 1.6 nm. 
The organic material was out-gassed in vacuum for sev- 
eral days at ~ 190 °C in order to remove residual rubrene 
peroxide [T3] , Well cleaned Si(100) substrates with a 
thermal oxide layer (thickness ~ 20 nm) were transferred 
into the growth chamber and kept at temperatures above 
400 °C for several hours. Before deposition of the organic 
material we determined the thickness of the silicon ox- 
ide and the precise angle of incidence by ellipsometry. 
The rubrene films were deposited at room temperature 
via sublimation from a Knudsen cell at 210 °C. A typi- 
cal deposition rate of 0.2nm/min and a pressure of less 
than 1 x 10 -9 mbar during growth ensured reproducible 
results. 

The experimental data were analyzed with the com- 
mercial WVASE software package [14]. Since rubrene 
grows in amorphous films, we used an isotropic model 
consisting of three layers, i.e. rubrene - silicon dioxide 
- silicon, with the optical constants of silicon and sili- 
con dioxide taken from Ref. [15] . First, we determined 
the thickness of the rubrene film by employing a Cauchy 
model in the spectral range below 2eV, where the ab- 
sorption of the organic material is negligible. Then, we 
performed a point-by-point fit which yields the complex 
dielectric constant e\+ie2 for each energy, see Fig.[2ja,b). 
To check the Kramers-Kronig consistency of the optical 
constants we also analyzed the experimental data with 
a general oscillator model. In this approach the imag- 
inary part of the dielectric function £2 is described by 
a sum of Gaussian peaks and the real part S\ is de- 
rived by the Kramers-Kronig relation. Generally, we 
found excellent agreement of the results obtained either 
by point-by-point fits or the general oscillator model. 
To check the consistency of the results obtained by UV- 
Vis spectroscopy and spectroscopic ellipsometry we also 
measured transmission spectra of rubrene on glass, see 
Fig. [2j;. We found that these data can be described very 
well using the optical constants of rubrene on silicon ob- 
tained by ellipsometry (and those of glass measured sep- 
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FIG. 2: (Color online) Complex dielectric function of a 33 nm 
thick rubrene film on Si02- (a) Real part e\ and (b) Imaginary 
part £2 • The optical properties were obtained in situ by spec- 
troscopic ellipsometry and analyzed using a point-by-point fit. 
(c) Light transmission through a 34 nm thick rubrene film on 
glass measured in air (solid black). A model, which is based 
on the optical constants of rubrene shown above (a, b) and 
the optical constants of glass obtained separately by variable 
angle spectroscopic ellipsometry, can be used to describe the 
data (dashed red). 



arately with variable angle ellipsometry). The minor dif- 
ferences between the data and the model might be due 
to oxidation of some rubrene molecules in air [16] [17] . 
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FIG. 3: (Color online) Displaced harmonic oscillator model. 
Because of different electron configurations in the So (HOMO) 
and the Si (LUMO) electronic excitations result in spatial 
deformations of the molecule. The corresponding change of 
the conflgurational coordinate Q gives rise to the observed 
vibronic excitations. 



III. RESULTS 



is virtually constant towards high energies and £2 is then 
proportional to a/E. Fig. [4^i shows the normalized data 
for rubrene in acetone solution and thin films together 
with the corresponding least-square fits. 

With the sub-band energies given by Eq. , their rel- 
ative intensities I n according to the Poisson progression 
in Eq. ([2|, and a Gaussian for each sub-band, the imag- 
inary part of the dielectric response can be expressed as 



62(E) 



5 



In 



exp 



(E - E n ) 
2al 



(3) 



where the broadening parameters a n are related to 
the respective widths of each sub-band by FWHM n = 
a/8 In 2 a n . Since several internal modes might contribute 
to each sub-band, the linewidth are not necessarily con- 
stant. In the fitting algorithm all parameters except the 
width of the last two peaks (fixed at a 4 = 0.08 eV and 
(J5 = 0.1 eV) have been varied. The Huang-Rhys factor 
S is then deduced from the intensities of the first four 
transition peaks using Eq. ([2|, see Fig. [27 . The re- 
sults of the fitting procedure for the vibronic progression 
are summarized in Tab. [U 





Eq-o 


E v 


S 


Solution 


2.359 


0.165 


0.986 


Thin Film 


2.335 


0.169 


0.974 



A. Analysis of the optical spectra 

The imaginary part of the dielectric function (Fig. [2]o) 
shows a pronounced vibronic progression with five clearly 
distinguishable sub-bands. The displaced harmonic oscil- 
lator model - a simplified scheme describing the dipole- 
allowed transitions between the electronic ground state 
(So) an d the first excited state (SV) including spatial de- 
formations of the molecule (Fig. [3| - has been used to 
model the data. In this approximation the excitation 
energies E n are given by 

E n = E _ + nE v , (1) 

where £?o-o is the energy of the Soo —> $10 transition, E v 
the vibrational energy and n the vibrational quantum 
number of the excited state. The intensities I n of the 
vibronic sub-bands contribute to the imaginary part of 
the dielectric function according to a Poisson progression 

on 

I n = I^-e~ s , (2) 
n\ 

where S is the electron-phonon coupling constant, also 
known as the Huang-Rhys factor. This coupling param- 
eter essentially determines the shape of the absorption 
spectrum, i.e. the relative intensities I n of the vibronic 
sub-bands, which are given by the area under each peak. 

In order to compare the data taken in solution with 
the thin film spectra the absorption coefficient a of the 
solution has been converted to the dielectric function us- 
ing the common approximation that the refractive index 



TABLE I: Parameters derived from the experimental data us- 
ing the model of the displaced harmonic oscillator for rubrene 
in acetone solution and thin films: transition energy between 
the vibrational ground states Eo-o (eV), vibrational energy 
E v (eV) and Huang-Rhys factor S. 



B. DFT calculations 

In order to study the effect of possible conformational 
changes of the molecule and to compare the experimental 
results with theory we have computed the excitation en- 
ergies and the deformation in the lowest excited state of 
rubrene using the turbomole 5.7 program package [18]. 
We have chosen the B3LYP hybrid functional since it 
gives excellent molecular geometries of organic molecules 
together with rather precise transition energies resulting 
from a compensation between the systematic deviations 
arising from the underlying gradient-corrected density 
functional and the Hartree-Fock contribution [19 . For 
the electronic orbitals, we use a triple-^ valence polar- 
ized (TZVP) basis [20] throughout this article. 

In order to cover the twisted and the planar isomers 
of rubrene, we have optimized both of them in their re- 
spective point groups, i.e. D2 for the twisted isomer, and 
C2h for the planar one, see Fig. [5] Our simulation re- 
produces the known fact that the twisted isomer is more 
stable [10], with an energetic difference of 0.172 eV at 
the B3LYP/TZVP level. 
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FIG. 4: (Color online) (a) Normalized dielectric function £2 
of a 33 nm thick rubrene film (solid black) and of a dilute 
solution of rubrene in acetone (solid dark blue) with a fit 
to the data (dashed, red and light blue) based on Eq. (|3j). 
The absorption spectrum shows an redshift which is caused 
primarily by the refractive index of the environment ('solvent 
shift'), (b) Least-square fit to the intensities I n of the vibronic 
progression using an interpolation of Eq. ^ which yields the 
Huang- Rhys factors S for the thin film and the monomer, (c) 
The width (FWHM) of the transitions peaks. In case of the 
solution it shows a monotonic increase, which indicates that 
several vibrational modes contribute to each sub-band. 



In solution the charge distribution of the solute 
molecule polarizes the medium surrounding it, result- 
ing in a further stabilization energy. This phenomenon 
can be handled with the conductor-like screening model 
(COSMO), where the solute is treated as lying inside a 
cavity surrounded by the polarizable medium [21]. In 
this approach, the response of the dielectric medium 
is replaced by an ideal conductor, generating screen- 
ing charges on the surface of the cavity. In a second 
step these screening charges are scaled by a factor / 
depending on the dielectric constant e of the solvent, 
f(e) — (e — l)/(e + |). We have applied the respective 
routine as implemented in the turbomole 5.7 program 
package to a geometry optimization of each isomer in a 
medium with a dielectric constant of epsilon=20.7 cor- 
responding to acetone. At the B3LYP/TZVP level we 
found stabilization energies of 0.312 eV for the planar 
isomer, and 0.303 eV for the twisted isomer. This is re- 
ducing the energetic difference between the two isomers 
by a small amount, but the twisted isomer still remains 
more stable by 0.163 eV. 

Even though the twisted isomer is more stable both 
as a free molecule and in solution, in the crystalline 
phase the planar isomer can be stabilized by a beneficial 
geometric arrangement between neighboring molecules, 
resulting in a large cohesive energy so that the ener- 
getic cost of the planarization is overcompensated by at- 
tractive inter-molecular interactions. In the amorphous 
phase of solid rubrene it is expected that the twisted ge- 
ometry is conserved because the random positions and 
orientations of the neighboring molecules cannot stabi- 
lize an unfavorable isomer. 

The different energies of both isomers raise the ques- 
tion if the planar isomer is stable, or if it can sponta- 
neously transform into the more favorable twisted ge- 
ometry. Using B3LYP and basis sets in the size range 
between 6-31 G(d) and TZVP we found that the lowest 
symmetry-breaking A u mode in the Cih point group re- 
ceived a small imaginary frequency. For the larger TZVP 
basis this apparent instability was found to depend on the 
integration grid: A calculation using the standard m3 
grid resulted in an imaginary frequency of zl3 cm -1 , but 
based on a denser m4 grid, this value decreased towards 
ib cm -1 . Using a tighter convergence of 10 -8 Hartree, 
the energetic minimum along this imaginary mode oc- 
curs at 12 cm -1 = 1.5 meV below the /^-symmetric 
isomer. However, independently of basis size and integra- 
tion grid, for larger deformations along this vibrational 
eigenvector, the shape of the potential returns to an es- 
sentially parabolic minimum, with a tiny distortion for 
small deformations. Hence, this ultra-soft or even imagi- 
nary mode does not constitute a viable barrierless path- 
way from the planar isomer to the twisted configuration. 
For the twisted isomer we found positive frequencies for 
all vibrational modes, the lowest being a breathing mode 
at 26 cm -1 . 

In the twisted isomer the repulsion between the phenyl 
side groups achieves a rather larger distance between 



them as opposed to the planar isomer, where a relatively 
small distance between the phenyl groups is enforced by 
the rigidity of the bond connecting each side group to the 
tetracene core and by the more restrictive point group 
C2h- Therefore, the twisted isomer gains a substantial 
amount of energy through a reduction of these repulsive 
interactions, allowing eventually to invest a part of this 
energy into the unfavorable twist of the tetracene back- 
bone. The resulting angle between the two central rings 
of the tetracene core is 22.8° in the electronic ground 
state, and the angle between the two final rings 42.0°, 
compare Fig. (5}d. In the relaxed excited state these an- 
gles increase to 26.7° and 43.6°, respectively. 

For both isomers the relaxed excited geometries are 
computed with time-dependent density functional the- 
ory (TD-DFT) at the B3LYP/TZVP level, conserving 
in each case the point group of the electronic ground 
state. The resulting deformation pattern is projected 
onto the symmetry-conserving breathing modes of the 
respective isomer in its ground state, defining in turn a 
Huang- Rhys Sj factor for each vibration Tiojj . In contrast 
to the lowest symmetry-breaking vibration of the pla- 
nar isomer the frequencies of the symmetry-conserving 
modes depend only weakly on variational basis set and 
integration grid, with the largest influence for the low- 
est breathing mode found at 19 cm" 1 (B3LYP/TZVP, 
ra3) or 22 cm" 1 (B3LYP/TZVP, ra4). Correspondingly, 
we have determined similar values for the Huang-Rhys 
factors obtained at these levels. Applying the projection 
scheme to the deformation pattern of a cationic molecule, 
we have reproduced the key features of previous calcu- 
lations addressing the Huang-Rhys factors in the planar 
isomer [22]. Of course, when compared to the cationic ge- 
ometry, a neutral excitation of a planar rubrene molecule 
produces a somewhat different set of Huang- Rhys factors. 

By averaging over vibrations in the range between 
900 and 1800 cm -1 we define an effective mode with a 
Huang- Rhys factor S e fi = J2j &j •> a reorganization energy 
^eff = J2j SjTioUj, and a mode energy of Tiuo e R = \ e fi/S e fi. 
The harmonic approximation to the potential energy sur- 
face results in slightly too large vibrational frequencies. 
In order to compensate for this deviation we apply a scal- 
ing factor of 0.973 adequate for the B3LYP functional to 
the vibrational frequencies [23j [24] , conserving however 
the reorganization energies Xj assigned to each mode, so 
that the Huang-Rhys factors are increased by a factor 
of 1.028 = 1/0.973. The results of this computational 
procedure are reported in Tab. [H] We use two kinds of 
relaxed excited geometries: First, the result of the TD- 
DFT geometry optimization, and second, a modified ex- 
cited geometry where changes of the twisting angle of 
each phenyl group around the bond connecting it to the 
tetracene core are eliminated. From similar investiga- 
tions on other molecules containing phenyl groups we ex- 
pect that this procedure minimizes artifacts arising from 
finite angle effects related to the relatively large changes 
of these angles in the relaxed excited geometry [25] . 

When basing our projection directly onto the TD-DFT 
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FIG. 5: (Color online) (a) Geometry of the most stable 
rubrene isomer with a twisted tetracene backbone, and (b) 
planar rubrene isomer, which resembles the geometry in the 
crystalline phase. 



geometry for the excited state, the sum over the reorgani- 
zation energies assigned to each mode exceeds the avail- 
able reorganization energies and A^ e ^ deduced from 
the potential energy surfaces of ground and excited state, 
respectively. Using instead the modified excited state 
geometry where the orientation of the phenyl groups 
is adjusted to the ground state geometry, this artifact 
disappears, indicating that the non-orthogonality prob- 
lems arising from the modified orientation of the phenyl 
rings are eliminated. Moreover, in the second projection 
scheme the frequency assigned to the effective mode is in 
better agreement with the observations. 





twisted 


planar 


i^abs 


2.173 


2.279 


^0-0 


1.965 


2.085 


X (9) 


0.204 


0.195 


A (e) 


0.208 


0.194 


huj eS (TD-DFT) 


0.155 


0.150 


S eff (TD-DFT) 


1.216 


1.358 


hcueff (no twisting) 


0.162 


0.158 


S e ff (no twisting) 


0.985 


1.025 



TABLE II: Vertical transition energy _E a bs (eV) in the ground 
state geometry, transition energy between lowest vibrational 
levels Eo-o (eV), reorganization energies A^ (eV) and A^ 
(eV) on ground and excited state potential energy surfaces, 
respectively, effective mode Tiuj e R (eV), and effective Huang- 
Rhys factor for the twisted and planar isomers of rubrene, 
before (TD-DFT) and after (no twisting) eliminating the in- 
fluence of a modified orientation of the phenyl side groups 
in the relaxed excited geometry. All B3LYP vibrational 
frequencies are scaled down by a factor of 0.973, and the 
Huang- Rhys factors S are scaled up correspondingly by a fac- 
tor of 1/0.973 = 1.028, so that the reorganization energies 
are conserved. All calculations have been performed at the 
B3LYP/TZVP level. 
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IV. DISCUSSION 

The spectra of rubrene in solution and rubrene thin 
films show a pronounced vibronic progression that we 
have analyzed using the displaced harmonic oscillator 
model. For evaluation of the results the parameters com- 
piled in Tab. |T] should be compared with each other and 
with previous investigations: We find that the transi- 
tion energy i^o-o = 2.359 eV of rubrene in acetone so- 
lution deviates less than 2 meV from a value reported 
earlier [26]. Based on the systematic study of solvent 
shifts presented in that work one can estimate that the 
HOMO-LUMO transition of rubrene dissolved in acetone 
is red-shifted by about 0.100 eV with respect to a free 
molecule. Consequently, the lowest vibronic sub-band in 
the free molecule would be found at 2.459 eV. This illus- 
trates that the transition energy £?o-o = 2.335 eV mea- 
sured for rubrene in thin films shows a significant redshift 
which is caused by the environment. Since the Huang- 
Rhys factors for rubrene monomers dissolved in acetone 
and amorphous thin films show only minor differences, 
i.e. S so \ = 0.986 and 5mm = 0.974, we conclude that 
the exciton-phonon coupling is not markedly influenced 
by aggregation of the molecules - most likely because the 
investigated thin films are amorphous and the phenyl side 
groups keep those neighboring tetracene backbones apart 
that are mainly involved in the 7r — > 7r* transition. 

Furthermore, we observe a continuous increase of the 
peak width in the vibronic progression of rubrene which 
indicates that more than one vibrational mode is in- 
volved. For the solution spectrum we find a simple mono- 
tonic dependency, whereas the thin film spectrum shows 
a different behavior with the width of the first peak being 
larger than the second one. One could speculate that this 
observation is related to the reported existence of a con- 
formational rubrene isomer with a slightly different elec- 
tronic transition energy [10] . In order to verify whether 
two different isomers contribute to the optical spectra 
one might compare the experimental data with the re- 
sults of the DFT calculations for the twisted and pla- 
nar molecule (Tab, [ill): Despite the excellent agreement 
of the calculated Huang-Rhys factors with our experi- 
mental values [28] their small difference (Stwist = 0.985, 
5 p ian = 1-025) cannot provide the required sensitivity to 
to distinguish both isomers. Our calculations, however, 
show that the transition energies of the twisted and pla- 
nar monomer differ significantly with A^o-o = 0.120 eV. 
Unlike spectra of the pure isomers a mixture of the two 
molecules therefore would have to show some unusual 
spectral features. 

Based on the solution spectrum of rubrene, which itself 
represents the twisted isomer, we also model the spec- 
trum of the planar isomer by blue shifting the spectrum 
by 0.12 eV. By adding both spectra with different weights 
and normalizing them we simulate a corresponding mix- 
ture of planar and twisted molecules, see Fig. [6] We 
find that the spectrum changes mainly in the low-energy 
range when the number of twisted molecules is increased 
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FIG. 6: (Color online) Simulation of the spectra for differ- 
ent mixing ratios of planar and twisted rubrene isomers. To 
compensate the solvent shift and to allow a comparison of the 
line shapes the thin film spectrum is shifted in energy such 
that the i^o-o energy agrees with the simulated spectra, (a) 
mainly twisted isomers (b) mainly planar isomers. 



from e.g. 10 to 30% (Fig. |6p), whereas in the other case 
increasing the number of planar molecules from to 10% 
affects mainly the high-energy range (Fig. . Overall 
the simulations indicate that a contribution of the second 
rubrene isomer is relative small, i.e. below 10%. Because 
the twisted isomer is energetically favorable we conclude 
that the thin film spectrum is dominated by the twisted 
molecule and a fraction of no more than 10% of planar 
molecules may contribute to the spectrum. The differ- 
ences between the solution and the thin film spectrum 
in the high-energy range, however, cannot be explained 
by contributions of a small fraction of planar molecules 
alone. Additional coupling in the thin film, like exciton 
transfer, might be the origin of this deviation. 



V. SUMMARY AND CONCLUSIONS 

The optical spectra of rubrene thin films and 
monomers in solution have been investigated experimen- 
tally by spectroscopic ellipsometry and UV-Vis spec- 
troscopy. The analysis of the vibronic progression using 
the displaced harmonic oscillator model has revealed the 
subtle differences between the monomer and the thin film 
spectrum. Additionally, time-dependent DFT calcula- 
tions have been performed for the planar and the twisted 
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isomer in order to get more insight into possible confor- 
mational changes of the molecule as they are known from 
the literature. The theoretical results, notably the excita- 
tion energies and Huang-Rhys parameters for the planar 
and twisted molecule, are in excellent agreement with the 
experimental data. Since the two rubrene isomers show 
different transition energies, the optical spectra may re- 
veal their existence. In fact a simulation, which studies 
the effect of mixing the two species, indicates that one 
isomer - most likely the more stable twisted isomer - 
dominates the thin film spectrum with a percentage of 
more than 90%. 
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